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^sj , The effect of hole doping on the Tight-Binding {TB) model of the Cu-O planes in the La2Cu04, 

constructed in previous works is investigated here. Firstly, it is pointed out that the model employed 
■ constitutes a generalization of the Hubbard one for the same system. Thus, the former predictions of 

CNj _ the insulator gap, antiferromagnetic (AF) character and the existence of a paramagnetic-pseudogap 

(PPG) state at half-filling, become natural ones to be expected from this more general picture. 
^ ■ The effect of hole doping on the antiferromagnetic-insulator state (AFI) and the paramagnetic- 

h^- I pseudogap one at half-filling, is investigated here at T = OK. The results predict the occurrence 

, of a quantum phase transition (QPT) from the antiferromagnetic- insulator state at low doping 

\^ ■ to a paramagnetic-metallic state (PM) at higher hole densities. Therefore, a clear description of 

^S) ' the hidden QPT laying beneath the "dome" in high critical temperature (HTc) superconducting 

materials is found. At low doping, the systems prefers the AFI state, and at the critical value of 
the doping density Sc = 0.2, the energy of a metallic state starts becoming lower. The evolution 
^ ■ with small doping values of the band spectrum of the AFI state, shows that the holes tend to 

become localized at the middle of the sides of the reduced Brillouin zone (BZ). Then, when 5 
passes through the critical value, the holes of the AFI state move to become situated at the corners 
^ ■ of the same reduced BZ, showing in this way a structural change at the phase transition point. 

Thus, the PAI state which appears at the transition point acquires the same behavior with respect 

e*"^ . to the localization of holes as the pseudogap state. In the small doping limit a clear difference 

between the degree of convergence of the iterative self-consistent solution is associated to an even 
or odd number of electrons. It suggests that the Kramers degeneration in combination with the 
spin-spatial entangled nature of the hole states, leads to a new kind of pair interaction between two 
holes. The binding energy value is estimated as a function of the screening. 
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I. INTRODUCTION 
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■ Since the discovery of the high critical temperature superconductivity in the called "cuprates" (made up of copper 

■ oxides) [H, the search about the nature of their behavior has been and it continues being, one of the most defiant 
. and investigated topics in Condensed Matter Physics. Among these compounds is the La2CuOi, that is a classic 
' example of the so-called strongly correlated electron systems (SCES). However, due to the inherent difficulties of 

ps) . the many body problem for this class of materials, numerous properties continue without being clearing up, in spite 
T-H ■ of the efforts done by the investigators in this subject. The La2Cu04 superconductor is characterized by a drastic 
^ ] change in its behavior depending on the density of electrons in the two-dimensional planes Cu02 2], and it happens 
• ^ ■ in general in the HTc superconductivity materials. All the cuprate superconductors have as a common element, the 
] existence in their crystal structure of Cu02 planes. 

The standard band calculations predict a character of paramagnetic-conductor to this material j31] , in drastic contrast 
with its experimentally observed, insulating and antiferromagnetic nature. These two properties are associated with 
the presence of strong correlations, and are not derivable starting from an independent particle scheme as the Hartree- 
Fock {HF) one, when the whole many electron crystal system is considered However, in Refs. [34^, on the base of 
a single band model solved in the HF approximation, it was possible to predict the mentioned SCES properties of this 
material. At first sight, it might seem that to obtain the AF and insulating properties of La2CuOi from a HF scheme, 
could be a contradictory outcome, since those properties are thought as essentially indicating the strong correlation 
nature of this material. However, those results become natural ones, when considering the following circumstances. In 
first place, the HF procedure was not applied in its quality of a First Principle method to solve the full Hamiltonian 
problem associated to the total electronic structure of La2CuOi. Alternatively, the procedure was employed for solving 
the simple model of the Cu02 layers built in Refs. 0-01 • It is well-known that the impossibility of describing strong 
correlation properties by employing the HF scheme, appears when the scheme is applied to find the solution of the 
exact many electron problem defined by all the electrons and nuclei constituting a crystalline solid. For example, the 
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direct HF calculation of the band structure of La2CuOi, predicts an enormous electronic gap of nearly of 11 eV [7|. 
In addition, the model introduced in Refs. only considers the electrons that half- filled the single band crossing 

the Fermi level appearing in the early electronic band calculation Q. In fact, the proposed simple Hamiltonian can 
be viewed as a preliminary stage within the definition of one among the variety of Hubbard's Hamiltonians. One, 
in which the last step of imposing the nearest neighbor approximations leading to the Hubbard theories had not 
been implemented. Therefore, the proposed model could be considered as an improvement of a Hubbard one, which 
retains the full Coulomb repulsion operator intact [s^. This last point can be understood by keeping in mind that 
its free Hamiltonian is basically a TB one in which the full Coulomb interaction operator between the electrons 
had been retained. Thus, after performing the nearest neighbor approximations, a typical Hubbard scheme arises. 
Henceforth, since the Hubbard model is recognized to convey strong correlation effects, and moreover, since it is 
also known that its mean field solutions also can exhibit such effects Q, the appearance of the insulator gap and 
the antiferromagnetic structure in Refs. becomes a reasonable conclusion. However, a new physical prediction 

coming from the analysis is the existence of pseudogap states, which are not following in the Hubbard approximation. 
These mentioned results emerged as a consequence of a combination of lattice symmetry breaking with a spin-space 
entanglement structure in the single particle solutions of the HF problem. The treatment gave a clear explanation 
about the nature of the antiferromagnetism and the insulating structure of this high temperature superconductor. 

The process through which the electronic structure evolves with the hole doping, from the AFI state at low doping 
to the superconductor state, and then to a normal metal phase at large hole concentrations, has been the subject of 
a strong interest in the literature. It is known that the magnitude of the gap in the normal state is of the same order 
that the superconductor gap Q. With the objective of studying the hole doping dependence of the band structure, 
the La2-xSrxCu04 (LSCO) system is for several reasons appropriate among the family of the HTc superconductors. 
In first place, the LSCO has a simple crystalline structure of Cu02 layers. Second, the concentration of holes in 
these planes can be controlled in a wide range and is uniquely determined by the concentration x of Sr Q . For this 
reason, one can develop samples of the material by varying doping continuously from the insulating state without 
doping (x = 0) up to the high doping limit {x — 0.3) in the same system. 

Numerous experiments made in copper oxide compounds, suggest the existence of a quantum phase transition at 
T — K, which lays inside the doping interval of the superconducting dome (this is colloquially referred as "laying 
beneath" the dome). It is believed that this QPT is the key to understand the high temperature superconductivity 
and also to explain the properties of the normal state in those materials fl^, [TTj . The present work presents results 
that predict the existence of this sort of QPT in the context of the model for the La2CuOi constructed in Refs. [4-^^ 
when the hole doping is incorporated. The HF system of equations is solved here by employing the same method 
used in Refs. [J-Q- present case, again the parameters of the model are determined of imposing the condition 

that the HF solution without crystalline symmetry breaking, reproduces the form of the single band crossing the 
Fermi level in the band calculations of the material 0. Then, we investigated the changes that occur in the band 
structure when the system is doped with holes. As in Refs. at exact half-filling, the HF solution predicts 

the existence of the AFI and PPG states. This previous result motivated the idea of investigating the hole doping 
consequences, after conceiving the possibility of predicting the quantum phase transition at T = . The same 
existence of this transition, is currently one of the most fundamental questions in high temperature superconductivity 
research [l^, [HI ■ A transition like this is now considered that should determine the properties of the normal state 
in the diverse regions of the phase diagram. The present investigation aim consists in studying the existence of this 
quantum phase transition beneath the superconducting dome. 

The Hartree-Fock description of the simple model for the Cu-O planes of La2CuOi including symmetry breaking 
effects and spin-spatial "entanglement" of the single particle states at half-filling, is here extended to consider the 
effects of hole doping. It is observed how the effect of doping is able to predict a variety of the most interesting 
properties of this material a. T = Q K . The evolution of the band spectrum of the antiferromagnetic-insulator and 
the paramagnetic-pseudogap states as functions of the hole doping parameter are determined, for a wide range of the 
hole concentration < a; < 0.3. Around the critical doping Xc = 0.2, the results show that for the AFI state, the 
band spectrum suffers a gradual change, in which the insulator gap diminishes to completely closing. Surprisingly, 
the same behavior occurs for the PPG state for which the pseudogap also collapses. For higher values than the 
mentioned critical doping one, the AFI and PPG states coalesce in a single HF solution. The magnetization of this 
high doping state vanishes as well as its gap, thus, it describes a paramagnetic-metallic phase. The results also reveal 
a drastic change in the Fermi surface, which goes from a hole-like Fermi surface centered at (tt, tt) for < x < 0.2 
into an electron-like one centered at (0, 0) for 0.2 < x < 0.3. These results evidence the existence of a quantum phase 
transition laying beneath the superconductor dome, in which an insulator state with antiferromagnetic correlations 
transits to a paramagnetic-metallic phase. The work also identifies a new mechanism of hole pairing which could give 
rise to the superconductivity. The effect results from a combination of the Kramers degeneration with the spin-spatial 
entanglement of the single particle HF states. The discussion helps to clarify the relation between: the pair of states 
AFI and PPG in the Physics of the strongly correlated electron systems with the superconductor state and the 
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quantum critical point, around which the phase transition occurs |10l - ll3l |. An evaluation of the hole pairs binding 
energy as a function of the dielectric constant is obtained, which predicts binding of the pairs at the experimentally 
measured values of the dielectric constant of La2Cu04 [ij]. However, the thermodynamic limit for this calculation 
had not been attained yet in this work. In addition, the parameters of the model should be optimized yet to match 
the measured parameters of the material as the insulator gap of 2 eV (isj , and the dielectric constant having a value 
of the order of 25 [l3| ■ Here, it is pointed out that this binding effect could be the acting mechanism determining the 
physical relevance of the doubly charged bosonic fields argued in Ref. [l^. 

The papers is organized as follows. In section 2 we review the one band model of the Cu-O planes in La2C'uOi 
introduced in Refs. the procedure for the determination of its parameters and the Hartree-Fock solution for the 

here considered situation: the study of the hole doping effects. In section 3 the results are presented for the evolution 
of the band spectrum of the PPG state in the range of doping < a; < 0.3. Next, the section 4 considers the same 
study for the AFI state. Section 5 presents the results for the changes of the Fermi surfaces as the doping is increased. 
The quantum phase transition properties are discussed in section 6. At the section 7 we exposes the identification 
of a possible hole pairing mechanisms which could give rise to the superconductivity and estimate the values of the 
energy for pair binding. Finally, the conclusion are presented in section 8. 

II. THE Cu02 MODEL INCLUDING HOLE DOPING 

In this section we will review the main ideas and elements defining the model for the La2CuOi introduced in 
Refs. In the first subsection, the structure and notation of the unrestricted Hartree-Fock scheme employed is 

described. Next, the model is presented. 

A. Fully unrestricted Hartree-Fock scheme 



The N electron system considered in the TB model discussed in Refs. [^-|6| was described by a fully unrestricted 
Slater determinant /„(a;i; si, xa?; sat) state constructed with single particle orbitals 4>ki(xi,Si) with i = l,...,iV, 
which arbitrarily depends on the spin variables at any point of the space. The index n represents the set of quantum 
numbers of the many electron system. In this subsection, as usual, generalized coordinates will be assumed to 
incorporate the spatial in common with the spin ones. For the electron case under study, the Slater determinant gets 
the expression 

f{xi;si,...,XN;sN) = 7=tE,,i,...,^„ <?^''?i (2^1' si)-<^'jjv(a;w>sjv), ^.-^^ 

Vj rji = fci, fcAf, 

where e''i' - '''" is the Levi-Civita tensor. The single particle HF orbitals satisfy a set of coupled integro-differential 
equations of the Pauli kind. This set is derived from the minimization of the HF energy functional of the system 
under the conditions of normalization of / and the normalization of all the single particle orbitals . For more 
details see Refs. In general, the Hamiltonian of usual electronic systems has a free term ^0(2:1) (a kinetic 

energy plus an interaction with an external field one) with an addition associated to the pair Coulomb interaction 
between the electrons in the form 

'H{xi,...,xn) = ^iioixi) + ]^^V{xi,Xj). (2) 

After performing the above referred minimization process by using the Lagrange multipliers scheme, a "fully unre- 
stricted" set of HF equations for the orbitals clik^ is obtained in the form 

)0^j(a; ,s )](l)r,{x,s)- 

ni s' 

(I)*^{X ,S )V{X,X )(j)riix ,S )](j)jj^{x,s) ^ er,(j)riix,s), (3) 

m s' 

where rj, 771 = ki, Un- In this completely unrestricted way, the HF set of equations was wrote by Dirac [l7l |. In this 
form, the orbitals are allowed to show arbitrary spin projection at any point of the space, and this determines the use 
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of the terms "fully unrestricted" for this set of equations. The total HF energy of the N electron system Ehf and 
the HF single particle energies have the forms 

Ehf = E^('?l'^o|'7) + 5E,,,,,i('7.'7i|^l'7i,'7) - 5E,,,,,i('?>'7i|^l'7>'?i>> 

Et) = (^/I'Hok) +E,,i(^?:»7il^l'7i,'7) -E^i(^>»7il^l'7,'7i)- 

The definitions of the basis states of the one electron band model and their scalar products appearing in the above 
equations are specified in Refs. 

1. The a and /3 spin constraints. 

The analysis in Refs. was initially motivated by the aim of examining the restrictions that could be introduced 
in the description of many particle systems, by the frequently employed assumption about that the HF single particle 
states should have a definite +1/2 {a) or —1/2 (/3) projection of their spin at all the spatial points. That is, to show 
the structure: 

^'=(^'^) = Uf(xW(.) type p. 

Whenever the spatial functions (f'ti^) identical, the HF evaluation is called a restricted one, and when they can 
be different it is described as an unrestricted one [l^. However, both of these cases are in fact restrictive ones for 
the allowed spin orientations of the HF orbitals. It was emphasized in references that these assumptions about 
the spin structure of the orbitals can be characterized as definite constraints which drastically limits the generality 
of the space of function in which the HF single particle orbitals are searched. It can be noted that the HF scheme 
introduced by Dirac [l3| does not include any restriction on the spin structure of the searched orbitals. 



B. The model for the CuOi planes 




FIG. 1: Band structure of the La2Cu04, calculated by Matheiss 1987. The single half-filled band suggests the 
consideration of the here investigated model, which is based in a set of interacting electrons subject to a crystalline 
TB potential, tightly binding them to the points of the lattice formed by the Cu atoms in the material. 

In figure [U a calculated band diagram of the La2Cu04 is illustrated. The band structure was obtained by the use 
of the LAPW (Linear Augmented Plane Waves) method ^S] . Let us describe below the formulation of the model for 
the Cu02 planes introduced in Refs. which is the starting point of the present work. It can be noted that the 

last band occupied by electrons is half-filled, a fact that indicates a metallic character for the material. The form of 
this band suggested the validity of an approximate TB description for the electrons populating it. The less bound 
kind of electron in the compound La2Cu04, is the unpaired one in the Cv?^ , which at variance with the in the 
plane, does not has its last shell [id) closed. These electrons, in a qualitative picture can be estimated to be the ones 
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constituting the mentioned unique band crossing the Fermi level in the band calculations done in the Ref. [3| . These 
electrons can be reasonably considered as tightly bound to the Cu atoms in the Cu02 planes, assumed the completion 
of the atoms in those planes. The last remarkSjSupport the idea of taking as the lattice defining the TB model, 
that one giving rise to the half-filled band in Ref. [3]. That is, the planar squared lattice of points coinciding with 
the sites of the Cu atoms in the Cu02 plane (see figure [2]). 

The presence of all the other electrons filling the rest of the bands, in common with all the neutralizing nuclear 
charges, were taken to play a double role in the model defined in Refs. [^-Q- In one sense, this system is assumed 
as a polarizable effective medium which screens the Coulomb repulsion by means of effective dielectric constant e. In 
second place, this system is assumed to generate a mean periodic Tight Binding potential W-y, which strongly confines 
the electrons of the single partially filled band to be close to the Cu atoms in the Cu02 planes. The model is assumed 
to be purely 2D, that is, the space coordinates are assumed to have only two planar components x = {xi,X2)- It 
is completed by assuming that the electrons in the partially filled band are also interacting with the potential Ft 
generated by a jellium charge density which neutralizes the net charge of the electrons of the model. The jellium 
charges were assumed as periodically concentrated and having a planar gaussian distribution with a radial width 
defined by the parameter b. In resume, the free Hamiltonian for the TB model for the electrons in the half-filled band 
has the form 

?io(x) = |^ + T4^^(x)+F,(x), (6) 
2m 

in which the periodic TB potential is satisfying the periodicity condition 



W^(x) = W^(x-HR), (7) 

and the jellium potential is defined by 

47reoe ^ J |x - y| 

in which the coordinate vectors of the planar Cu atoms are 

R — nxiP^xi +f^a;2^'6a:2 where nx-i^,nx2 G (9) 

where the unit vectors e^-j^ and e^.^ lay on directions defined by the vectors joining the nearest neighbors of the Cu 
atoms in the lattice in figure [2] It is known that the distance between a Cu atoms and its nearest neighbor is p = 3.82 
A [l^. Further, the Coulomb interaction between two electrons in the partially filled band of the model is assumed 
in the form 

1 

V=- r, (10) 

Aireoe |x - y| 

including a dielectric constant which is supposed to be generated by the polarization of the set of electrons filling the 
other bands and all the nuclei in the La2CuOi, through which the electrons of the partially filled band of the model 
are assumed to move. 

The model was started to be constructed with the initial idea of searching for HF single particle states of the 
electrons in the partially filled band, not being of an a or /? types, that is, not being separable in their spatial and 
spin dependence. In other words, allowing that in those HF orbitals the spin projection could depend on the spatial 
position (spin-space entanglement or non collinear spin structure). Materials showing AF structure were suspected 
to lead to HF solutions of this kind, a fact that also suggested the possibility of improving the understanding of the 
troublesome state of knowledge of the band structure of such materials. In the particular case of the La2Cu04, it is 
known that it is an antiferromagnetic compound. In the AF systems, normally the crystal symmetry of the substance 
is directly broken by the AF order. Therefore the description of the model was also designed to incorporate the 
possibility of the breaking of the crystalline translation invariance. The AF order has translation symmetry in each 
one of the two interpenetrating squared sublattices in which the planar Cu lattices can be decomposed. Thus, the 
HF single particle orbitals in the model were assumed to retain the Bloch functions character under the common 
reduced translation symmetry group of these sublattices. Then, the physical state should be equivalent under the 
translations leaving invariant the sublattices, but not under the one which transforms a sublattice into the other. This 
reduced symmetry is associated to a subgroup of the total translation group of the crystal and therefore its space 
representation should be more reduced in its number of states as classified by the momenta k in the reciprocal space. 
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FIG. 2: (a) The picture illustrates the Cu02 planes in the material, (b) Point lattice associated to those planes. In 
the search of the possibility of describing the AF properties of the material through the elimination of the symmetry 
restrictions (over the space of functions in which the single particle states were searched) it was helpful to 
decompose the full lattice in the two sublattices which are differentiated in the picture. 



The two sublattices of point depicted in figure [5J that will be described below by the indices r 
analytically defined in the form 



R('') = V2nipqi + V2n2pq2 + q' 



(r) 



1 or 2, are 



(11) 



if r = l 
pex, i/ r = 2, 



where qi and q2 are the two unit vectors defining the directions of the unit cell vector of the sublattices. 

1. Sublattice translations 

Then, the searched HF single particle wave functions should be eigenfunctions of the discrete translation operator 
Tj^(i), which transform a sublattice in itself 



(12) 



When the point lattice is infinite, the Brillouin cell formed by the set of all the momenta k indexing the eigenfunctions 
of Tp^{i) is the shadowed zone in figure [31 while the large square represents the set of momenta associated to the 
eigenfunctions of the group of translations in the original lattice. Since it is impossible to numerically treat the 
infinite lattices because it has a continuum of states, periodical boundary conditions were chosen in order to make 
finite the number of eigenfunctions in the Brillouin zone. The boundary conditions were imposed fixing the periodicity 
of the eigenfunctions (jj^^i in the boundaries defined by xi — —Lp and xi — Lp along the xi axis and by X2 — —Lp 
and X2— Lp along the X2 axis. The finite set of momenta vectors obeying the boundary conditions are 



e Z 



2 — ^xi i ^X2 ^ 2 ■ 



The number of electron eigenfunctions obeying the conditions is -^i that is, a half of the number of functions satisfying 
the same boundary periodicity conditions when the group is the group of translations in the full point lattice. 



2. Ttght-Btndmg single band basis 
The TB basis that was employed for describing the one band model had the form 



(s) Vexp(*k.RW)^o(x-RW), 



R('-) 



(13) 
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FIG. 3: The darker region shows the reduced Brillouin zone associated to the Bloch functions defined on the 

sublattices. The length of its sides is -y/Z^- 



V5o(x) 



exp(-^), a<p. 



(14) 



where N is the number of electrons in the model electron gas filling part of the band crossing the Fermi level and 
is the spin projection operator in the z axis, which is assumed to be orthogonal to the Cu02 planes ; cr^ = — 1 , 1 
are its eigenvalues; r = 1 and 2, is the above defined index for each of the two sublattices. Note that in the 
small overlapping approximation between nearest neighborhood points (which should be in different sublattices) , the 
orthogonality between the members of this basis is only lost between functions belonging to different sublattices having 
the same spin polarization and momentum. The orthogonality between distinct elements corresponding to the same 
sublattice is valid by construction. 

The chosen Wannier orbital of the model ipo{^ — R'-'^-') were fixed for simplicity. Their gaussian form implies that it 
was supposed that the net TB potential in the close neighborhood of each Cu atom had been assumed to be a simple 
2D harmonic potential. This simplifying assumption was adopted following the idea that the main forces in defining 
the AF and insulator properties of the Cu02 planes were in fact determined by the spontaneous breaking of the 
crystalline symmetry plus the spin entangled structure of the HF electron orbitals. Thus, simplifying considerations 
led to fix the explicit form of equation (IT^ for the starting basis states of the model. However, taking into account 
improved definitions for the Wannier's orbitals defining the single band model, might be of help for the description of 
other effects, by example the magnetic anisotropy of the AF order in La2CuOi. For this purpose, the model could 
be generalized by considering atomic ZD representation, employing Wannier's orbitals like the d ones of the Cu and 
the inclusion of spin-orbit interactions. This study is expected to be considered elsewhere. 



C. Hartree-Fock solution with hole doping 



In this subsection we present the matrix problem to which the HF set of equations of the model was reduced in 
Refs. after its projection on the defined basis. In the present work, these equations will be solved for the more 

general situation in which the system is doped with holes. In the coming subsection it will be reviewed how the 
model reproduces the main characteristics of the dispersion of the single band crossing the Fermi level in Ref. Q 
(predicting a metallic and paramagnetic state) when full crystalline symmetry and the HF orbitals of the type a 
and j3 are assumed. Afterwards, the next sections will illustrate the consequences of the elimination of the symmetry 
restrictions in the space of functions in which the HF energy functional is minimized. It was performed in a similar 
way as it was done in Refs. but for the important situation in which the systems is doped with holes. At very 

small doping values, the solutions are the ones obtained in Refs. for the exact half-filling situation: the ground 

state is an antiferromagnetic-insulator and the excited phase corresponds to the paramagnetic-pseudogap state. When 
the hole doping is augmented the results indicate the appearance of a quantum phase transition at the doping value 
X = 0.2, in which the AFI and PPG states both coalesce in one single metallic state for higher doping values. 
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As in Refs. [3-^ the searched HF single particle states are expressed as a linear combination of the before defined 
basis functions in the form 



(15) 



where I is the index of the quantum numbers required to uniquely define the HF single particle states. After 
substituting the above expression for the orbitals in the HF equations and taking their scalar product with an 
arbitrary state, the equivalent self-consistent matrix problem can be written in the form 



[El + x(G^° - GJT* - Fk)] • S"^-' = e,(k)/k • S"^-', 
with the definitions for the constants 

9 9 

me a 



(16) 



(17) 



e;(k) 



ma 



■ei(k), 



(18) 



which are dimensionless, as also are all the implicit parameters in the definitions of all the entering matrices 



El 



'^k,(t,r,a^,cr^) 



4x4 



G 



c 

k, (t,r,a2 z ) 



4x4 



Gl 



G 



Fk = Fi 



k,(i,r,a2,o'5 



I4x4 



-^fc ||-^k,(i,r,as,crs) 114^4 ■ 

They are respectively associated to the kinetic term plus the periodic TB potential, the Coulomb direct and exchange 
potential terms, the potential generated by the compensating jellium charge density and the overlapping matrix 
between the basis functions. The formulae for the matrix elements are shown in the Appendices of the Refs. @, @]. 
In the employed representation the normalization condition of the HF single particle states and the formula for the 
HF energy take the forms 



1 



HF 
k,/ 



(£F-e,(k))[ei(k) - ^B^''^*.{G^ - GD.B 



k,ii 



(19) 
(20) 



The HF matrix system of equations (|16p is a non linear one in the variables i?^ 'i^ : which are the four components 
of the vector B^'^ for each value of k. The B constants can be interpreted as determining the probability amplitudes 
to find the electron in the states (k, I), of the sublattice r, and spin along the z axis. Note that for each k value, four 
values of the quantum numbers of the HF single particle states I ~ 1,2,3, A should be obtained, or what is equivalent, 
four bands in the reduced BZ. From the equations (|16p it can be observed that the full matrix representation of the 
Fock operator is block diagonal in the momenta indices k, which is a direct consequence of the symmetry with respect 
to the reduced translation invariance of the system. 
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D. Maximal symmetry solutions: fixing the model parameters 

In first place the HF problem was solved by assuming that the HF single particle states satisfied the translation 
invariance in the whole lattice and also were showing a or /3 spin dependence. After this assumptions, as it was 
suspected, the HF solution of the model produced the single particle spectrum depicted in the figure HI in a similar 
way as in Refs. before we will adjust the free parameters e: the dielectric constant of the effective medium; 

a: the radial distance for which the gaussian Wannier's orbitals are appreciably different from zero; 7: the hopping 
probability between nearest neighbor sites fixed by the effective medium and b: radial distances inside which the 
jellium charges are concentrated. 

The Bloch Tight-Binding single band basis for this problem had the form 

(x, s) = J ^u^^ (s) exp(«Q • R)¥'o(x - R), (21) 

R 

where the appearing momenta Q are given by 

2 — 1 ^ "2" 

These functions define the Bloch states forming the basis of the maximal translation group with periodic boundary 
conditions in the same region defined before. In addition N = L x L, and R are the number of cells in the total lattice 
of points and the vectors defining the positions of the lattice. Let the HF single particle states expressed in the form 

0Q,(x,s)=^i??;VQ^(x,s), (22) 

in the just defined basis. In this case, the equivalent HF matrix problem obtained after substituting these functions 
in the set of HF equations, is now of dimension two. That is, there are two components I = 1,2 for each momentum 
Q. In this way, in an analogous form as in (I16p is possible to obtain the HF matrix problem in the form 

[E'q + X(G§ - - Fq)] ■ B^-^ ^ eKQ)/Q • \ (23) 

which is a system of non linear matrix equations to be solved by iterations. 

For starting the iterative process an initial paramagnetic state is employed as an ansatz. The figure |4] shows 
the doubly degenerated metallic and paramagnetic band which is obtained at half-filling conditions. That is, with 
iV = 20 X 20 electrons. The chosen parameters were: e = 12.5, which is a typical value in semiconductor systems. In 
what follows we will explain the reasons for selecting the values d — 0.25, b = 0.05 and 7 = 0.03, always following 
the criterium of obtaining a band width of 3.8 eV as it is suggested by the band diagram in figure [TJ 



620!l2o(eV) 




FIG. 4: The doubly degenerated paramagnetic and metallic band arising from the HF solution with full translation 

symmetry and a and /3 spin structure of the single particle states. 

Observe the coincidence in form between the calculated band and the conduction one appearing in figure [TJ In both 
of them, the Fermi level is crossed by the dispersion curve at the mid point between the top and the bottom of the 
bands, in the direction T-X, while in the direction laying at 45° respect to T-X the Fermi level becomes tangent to 
the dispersion curve in its maximum. 
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III. EVOLUTION OF THE PARAMAGNETIC-PSEUDOGAP BAND SPECTRUM 

In this section we will present the results of solving the HF matrix set of equations for a variable hole doping in the 
case of the excited solution which evolves from the pseudogap state obtained in Refs. [^-Q for the exact half-filling 
situation. It is important to recall that this state emerges as a solution after only eliminating the restriction on the 
single particle states of being a Bloch function in the full lattice formed by all the planar Cu atoms within the Cu02 
planes. The constraints of being single particle states of a or /3 types were yet maintained. This process allows to 
obtaining of a doubly degenerated HF solution which shows a pseudogap resembling the one observed in the normal 
state of the HTc superconductor materials. 

Figure [5] shows the evolution with increasing hole concentration of the bands of the pseudogap state, obtained for a 
point lattice of 20 x 20 points. All the graphics are plotted in the Brillouin zone of the sublattices, that is the darker 
zone in the figure [H Note the existence of a pseudogap which attains its maximum of the order ^ 80 meV. The 
parameters fixed in the previous section were employed for this evaluation. The maximal value for the gap appears 
at the mid points of the sides of the Brillouin zone of the sublattices and furnishes an estimate of T ~ 1000 K for the 
temperature at which the pseudogap starts to be observed in the experiments [13, HO, HII • 



e20x2o(eV), e20x2o(eV), x=Q.05 




FIG. 5: The evolution of the band spectrum of the PPG as the hole concentration grows. The Fermi level lays in 
the zero of the graphics. It can be observed how the pseudogap starts closing at the corners of the Brillouin zone as 

the doping increases. 

The evaluated result for the pseudogap depends from the effective dielectric constant e and the set of parameters 
which were fixed in order to reproduce the single band crossing the Fermi level in the band calculation done by 
Matheiss Then, the value obtained here should not be taken as precise prediction for the temperature of the 
observed pseudogap. However, the ARPES experimental results for the doped La2Cu04 indicates a pseudogap 
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temperature in the region 100-200 meV, that is value close to 1000 K. Thus, the PPG Hartree-Fock sohition offers a 
reasonable estimate of the pseudogap temperature T* [l^, [H, [13, HH . 

One interesting results which can be observed from figure [5] is that by augmenting the hole concentration, the 
so-called "pockets" [23 - |23 | in the corners of the Brillouin zone are formed. In addition it can observed that the 
pseudogap starts to diminish first in the corners until it fully collapses at the critical doping Xc = 0.2. After this 
critical concentration of holes, the state starts behaving as a paramagnetic metal. 



IV. EVOLUTION OF THE ANTIFERROMAGNETIC-INSULATOR BAND SPECTRUM 

In this section we will consider the solution of the HF system of equations (fTH]) by again using the successive 
iterations method, but in this case, starting from an antiferromagnetic initial ansatz for the Slater determinant state. 
The results presented below, were found for the values of the set of free parameters: e = 12.5, a — 0.25, 7 — 0.03 and 
b = 0.05, which were also employed in the previous section. 

The figure [S] (a) shows the band obtained for the point lattice of 20 x 20 points. The bands are clearly associated 
to an insulator system, which resulted as the ground state with respect to the pseudogap phase. 



e20s2o(eV), x=0 e20x2o(eV), x=a05 




FIG. 6: Evolution of the band spectrum API as the hole doping increases. The Fermi level is laying in the zero of 
the plots. The hole states form the the so-called Fermi arcs at the middle of the sides of the reduced Brillouin cell 
and the gap diminishes until it is fully closed at the critical doping Xc = 0.2. 

It also shows a local magnetization. These results are essentially the same as the ones obtained in Refs. '3-^6'] but 
for slightly different values for the free parameters of the model. The next figures |6] (b) , (c), (d) show the evolution 
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with the increasing hole doping of the band spectrum of these states. It can be observed that the added holes start 
accumulating at the mid points of the reduced Brillouin zone. That is, in the more energetic electronic states, and 
as a consequence, more preferable by the holes. All the pictures are plotted in the reduced Brillouin zone of the 
sublattices. When the doping is further increased up to some critical point, the holes, which were sitting at the mid 
sides of the Brillouin zone, drastically pass to occupy the corners of the reduced Brillouin zone. This constitute a 
structural change of the band occurring at the critical value of the hole doping 6c = 0.2, in which the form of the 
bands becomes different before and after the states passes through this particular hole concentration. In other words a 
quantum critical phase transition occurs at this value of the hole do ping . Thus, the model predicts a phase transition 
which had being argued to exist beneath the superconductor dome [lOllTl| . 

The results state that at the moment in which the AFI and PPG become degenerated at the quantum critical 
point, the obtained HF solution of the problem becomes a unique metallic and paramagnetic phase. An important 
point to underline, is that, as it is shown in figure [71 around the critical doping Xc = 0.2, the results predict that the 
insulator gap of the AFI state diminishes until it completely collapses. The same occurs to the state PPG. This 
produces the metallic behavior which coincides with the one observed in the material in the high doping region. 

Decrease of the insulating gap Collapse of the PPG gap at the nodal reg ion (x^) 




FIG. 7: Diminishing of the gap with the increase of the hole concentration, a) For the antiferromagnetic-insulator 
spectrum, b) For the paramagnetic-pseudogap state. Note that in both cases the gaps are completely closed when 
the concentration of holes is Xc ^ 0.2, predicting a metallic behavior. 



V. EVOLUTION OF THE FERMI SURFACE 

Let us expose now the results for the modifications of the Fermi surface induced by the hole doping. The Fermi 
surface, as separating in the reciprocal space, the occupied from the empty orbitals (25j |. allows to define many physical 
properties of the materials. It is known that in the largely non understood normal state of the HTc superconductors, 
the Fermi surface becomes truncated in parts called Fermi arcs. One of the most important open questions in 
superconductivity theory is how the Fermi arcs and the superconductivity are mutually related [261 127| . We have 
investigated the Fermi surface predicted by the model and its dependence on the hole concentration for a wide 
range of concentrations < 2^ < 0.3, defining the evolution with doping from the AF-insulator region up to the 
paramagnetic-metallic one. 

The figure [S] shows the evolution of the Fermi surface as the hole doping grows. It can be observed that the model 
predicts that in the low doping region (x < 0.15) the Fermi surface is composed of the above mentioned Fermi arcs 
around the nodal region (tTjTt). In this direction is where the doped holes establish as it was mentioned before (the 
grey zone in figure [5]). It was possible to evidence that the length of these arcs grows in proportion with the doping, 
up to the attainment of the critical doping at Xc — 0.2, where a sort of hole pockets form in the corners of the reduced 
Brillouin zone. 

In accordance with our calculations shown in figure [U the Fermi surface for x — 0.3 becomes squared and has a 
large straight portion around the point (-1, The areas closed by the Fermi surfaces are 70, 80, 85 and 95% of the 
half Brillouin zone area and corresponds to a; = 0.3, 0.2, 0.15 and 0.05, respectively. These proportions are consistent 
with the Luttinger sum rule, if the electron density is 1 — .t = (70, 80, 85 and 95%, respectively). It can be concluded 
that our results have a close qualitative agreement with the ones shown in figure 9, that were obtained by means of 
the ARPES in Ref. [9|. 

Therefore, it follows that the considered model predicts that the Fermi surface of the LSGO undergoes a transition, 
which goes from a hole-like Fermi surface centered at (tt, tt) for < a; < 0.2 into an electron-like one centered at 
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FIG. 8: Evolution of the Fermi surface for increasing hole doping. The areas closed by the Fermi surfaces are 70, 
80, 85 and 95% of the half Brillouin zone area and corresponds to x = 0.3, 0.2, 0.15 and 0.05, respectively. These 
proportions are consistent with the Luttinger sum rule, if the electron density is 1 — x = (70, 80, 85 and 95%, 

respectively). 



(0, 0) for 0.2 < X < 0.3. It is believed that the drastic change of the Fermi surface can correspond with the fact that 
the sign of the Hall coefficient changes from positive to negative values around x = 0.25 in the LSCO [1^ [2^. In 
the framework of the investigated model, the mentioned change is associated to the occurrence of a quantum phase 
transition around the critical point Xc = 0.2 which lay beneath the superconductor dome. 



VI. THE QUANTUM PHASE TRANSITION 



One of the current questions in HTc superconductor theory is the occurrence of the superconductivity as determined 
by the existence of a quantum critical point laying within the hole doping interval at which the superconductivity 
occurs. Today is widely debated the question about what is the detailed connection between the critical point and 
the superconductivity effect. Experimental indications about the existence of a critical quantum point comes from 
the transport and thermodynamic measurements [30l - l33| . In the framework of the present study, the HF energy per 
particle for the AFI and PPG states were evaluated by varying the hole concentration in the range < a; < 0.3, by 
using the formula 

E^f - E 0(e>-e,(k))[£Kk) - ^B^''*.{G^ - Gl).B^'% 
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FIG. 9: Evolution of the Fermi surface for increasing hole concentration for the La2~xSrxCu04,, obtained by the 
ARPES experiment. Taken from: Ino A 1999. The dotted line in the first figure represent the so-called Fermi arc. 

That is a discontinuous Fermi surface. 

The figure 10 (a) show how the AFI state, which have the lowest energy at half-filling, evolves and becomes 
degenerate with the PPG state at the critical doping Xc = 0.2. At this point, as it can be seen in the figure 6, 
the results predict for the AFI state, that the insulator gap diminishes to be completely closed. The same behavior 
is shown by the PPG state for which the momentum dependent pseudogap collapses. The solution of the problem 
determines a metallic behavior which is the observed nature of La2GuOi for high hole concentrations. 

Transition AFA-PPG 
-0.4 



(a) 




FIG. 10: (a) Dependence of the HF energy per particle with the hole doping for the AFI and PPG states, (b) 
Generic phase diagram where it can be observed the quantum phase transition from the AFI state to a 

paramagnetic-conductor. 

The magnetization of the AFI state was also evaluated by its defining formula 

"^W = "^k.ilx, s)a{s, s')0kj(x, s'), (24) 

k,/ s,s' 

where 

cr(s,s') = cr^i(s, s')ej.j -|- cr^2 (s, s^e^.^ -f ^^(s, s')ez, (25) 

and , (7x2 and are the Pauli matrices. One important conclusion in this study is that in approaching the critical 
doping value the local magnetization of the AFI state tends to vanish (see figure [TT|) , a property that the PPG 
showed along its whole hole doping evolution from the initial half-filled state. 

VII. A NEW HOLE BINDING MECHANISM: SPIN-SPACE ENTANGLEMENT 

The evaluation procedure of the energy for the AFI state in the region of low doping showed a clear dependence 
on the odd or even character of the number of holes added. That is, it was needed a very much larger number of 
iterations for attaining convergence when the number of holes was odd. 
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0.00 0,06 0.10 0.16 0,20 

Hole concentrations 

FIG. 11: Reduction of the magnetic moment of the AFI state under the doping increasing. 



We estimate that a reason for such a behavior might be that adding an odd number of holes, could be expected 
to show a larger complexity of the HF single particle states, since the Kramers degeneracy should be implemented. 
That is, it could be expected that the double degeneracy implied by the time inversion invariance of the system, when 
the odd number of holes does not allow to fill an integer number of degenerate pairs of states, should make such single 
particle states to show a very much complex structure. 

A perhaps related circumstance is that the spin-spatial entangled structure of the HF single particle states leads 
to imagine a mechanism of hole pairing that could be acting in the considered system. The idea is that this entangled 
structure of the single particle states might produce that two, well separated in space, hole wave packets could show, 
each of them, their relatively complex spin-spatial entangled constitution. However, when the two holes are allowed 
to be close in space, it is possible that they could tend to compensate their spatially dependent spin and magnetic 
moment structures and show a lower energy than the other pair of well separated holes. In the case of the paramagnetic 
systems such an effect seem to be very much weaker, due to the simplicity of the spin structure. 

In order to check for this possibility we defined a quantitative measure of the binding between two holes in the 
antiferromagnetic state. This definition was based in the energy of the ground state of the system Eq at half-filling 
and the ground state energies after doped with one hole Ei and with two holes E2 in the form 

A_B = 62 - 2ei, 

where ei ~ Ei — Eq and 62 — E2 — Eq. Whenever, two holes minimize their energy by producing a bound state, 
then Ab becomes negative. When vanishes the holes may not form a bound state since then 62 = 2ei and it is 
expected that they behave as independent excitations. 

The binding energy of two holes As as a function of the dielectric constant e for point lattices of 16 x 16 sites 
(squares) and 20 x 20 sites (circles) in the case of the AFI state is shown in figure [12] Note that the bound states of 
the two holes form at low doping value as helped by the amount of dielectric screening of the Coulomb interaction. 
The figure [12] indicates that in general, when the dielectric constant e increases, the binding energy reduces until they 
form a bound state. The critical value for the appearance of binding for the holes slowly grow with the increasing of 
the size of the sample (the region in which periodic boundary conditions were imposed). 

The results indicate that although the thermodynamical limit is not yet attained for the evaluation of this quantity 
Ab, it seemingly exists. An interesting result is that in La2CuOi, it has experimentally determined that a supercon- 
ductor gap Asc is the order of the 10 meV @, which is of the similar magnitude of the shown values in figure [12] for 
the two holes binding energy. That is, the measured superconductor gap approximately coincides with the energies 
A B required to break the two bound holes. This outcome suggests the possibility that the HF solution is able to also 
convey a pairing effect determining the presence of HTc superconductivity and the constitution of the Cooper pairs. 

A possibly acting binding mechanism could be as follows. Firstly, recall that the HF single particle states show 
the more complicate than the standard spin-spatial entangled structure. This complex composition can be expect 
to provoke, if the screening is assumed to be strongly enough, a diminishing energy effect, in which the respective 
entangled magnetic moment structure of the combined two states could tend to compensate one to another. Such 
an effect should appear in the bound state Bethe-Salpeter equations for two holes. In the described picture, the 
bound state of two holes at low doping in the AFI phase will be formed in this picture thanks to the screening of the 
Coulomb interaction. The fact that the measured static dielectric constant of the La2CuOi is around the value 25, 
indicates that the amount of screening is high. 

In the more investigated models, like the t — J one, it is argued that the existence of each hole breaks four 
antiferromagnetic bounds 0, an effect that have an energy cost of the order of the coupling energy (0.1 eV). Therefore, 
at least in the low hole density limit, like the case under consideration, two holes minimize the energy to create them 
by sharing a common bound. In this way they minimize the number of broken antiferromagnetic bounds. Note that 
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FIG. 12: Binding energy As of two holes as a function of the dielectric constant e for point lattices of 16 x 16 (blue 

squares) and 20 x 20 (red circles). 

this picture seems to be compatible with the one previously exposed on the basis of the spin-spatial entangled nature 
of the hole states. It is also known that the size of the Cooper pairs of the material is a small quantity (around two 
lattice constants) 0. 

It had been also argued that the existence of preformed hole pairs could not be sufficient evidence for the appearance 
of the superconductivity. In the context of the present model, it has been argued that for large values the dielectric 
constant and the presence of spin-spatial entanglement, the pairs could be formed. Then, let us assume that their 
pair wave function can be formed by superposing products of single particle waves having momenta values close to the 
centers of the four reduced Brillouin zone sides in figure 6. Then, the amount of the momentum transfer associated 
to these functions in the Bethe-Salpeter equation could be expected to produce pair wave functions showing sizes of 
the order of few lattice cells. If such results to be the situation, then, the hole pairs could be expected to condense 
at r = if giving rise to a Bose condensate of Cooper pairs showing superconductivity. These possibilities will be 
investigated elsewhere. In ending this section, we also want to express that the results also suggests that the here 
argued binding effects between two added holes, could constitute a dynamical foundation of the physical meaning of 
the charge two boson excitations identified in Ref. [l^, as being relevant for the description of Hubbard models. 

VIII. CONCLUSIONS 

We have investigated the model for the La2CuOi defined in Refs. applied to the situation in which the material 
is doped with holes. The main elements of the model were reviewed by fixing its free parameters and solving the 
Hartree-Fock equations. The evolution with doping of the band spectrum of the antiferromagnetic-insulator state and 
the metallic-paramagnetic pseudogap one are presented for a wide range of hole concentrations < a; < 0.3. Around 
the critical doping density Xc = 0.2, the results show that for the AFI state, the band spectrum suffers a change in 
structure of the bands and the insulating band gap diminishes until its complete closing. A similar process also occurs 
for the PPG state where the pseudogap collapses at the same special hole density. Over this critical density the 
system only shows a single paramagnetic-metallic state. The magnetization of the AFI state becomes zero exactly at 
the critical point. The solutions evidence that the destruction of the antiferromagnetic order is produced by the fact 
that added holes tend to occupy the states which show a more intense AF order: that ones which are closer to the 
boundary of the reduced Brillouin zone. The results show a drastic change of the nature of the Fermi surface, which 
goes from a hole- like Fermi surface centered at (tt, tt) for < a; < 0.2 into an electron-like one centered at (0, 0) for 
0.2 < a; < 0.3. Henceforth, the whole discussion had shown that the investigated model predicts the existence of a 
quantum phase transition at critical doping value which is beneath the superconducting dome in La2CuOi. It was 
also identified a possibility for the pairing of holes which could give rise to the superconductivity and the bounding 
energy was estimated. The effect seems to be closely linked to the spin-space "entanglement" effect of the HF single 
particle states. Possible connections with the existence of bosonic doubly charged excitations investigated in Ref. 
were also pointed out. 

In ending, let us comment on some possible extensions of the present work. It seems convenient to consider a new 
parameter for to be fixed in the Hamiltonian of the model. It is related with the effective mass defining the kinetic 
term of the free part of Hamiltonian. This constant was assumed to be equal to the electron mass in free space, which, 
possibly, is a somewhat rude choice, since such electrons are assumed to move in the crystalline potential generated 
by all the other particles, filling the rest of the many bands. This new freedom in the parameters might be relevant 
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when the spin-orbit effects wih be taken into account in order to describe the magnetic anisotropy of the AF order 
in further discussions. It would be also of help in allowing to more precisely fix the value of the insulator gap of the 
AFI state which is known to be close to 2 eV. One important point in this respect is that the dielectric constant e has 
been measured for La2Cu04. Therefore, it will be needed to fix the free parameters in order to define the observed 
value of the gap. At this point it can be noticed that electron effective mass will be an influential value which can 
control the energy scale of the bands and could be phenomenologically fixed to define the 2 eV insulator gap. It can 
be observed from equations p7)) and (IT51) . 

It seems also of interest to perform calculations in order to attain the thermodynamic limit in the hole pairing effect, 
and also to attempt deriving this binding effect but in the framework of the Bethe-Salpeter equation for two holes 
moving in the medium. This would confirm the presence of preformed Cooper pairs in the AFI state of the model. 
Lastly, it will seem convenient to improve the formulation of the model by employing 3D Wannier states resembling 
the incomplete 3D shell of the Cu atoms within the Cu02 planes. Upon this, the incorporation of the spin-orbit 
interactions will be needed in order to further describe the magnetic anisotropy of the AF order in La2CuOi. 
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